Are Gamma-Ray Bursts a Standard Energy Reservoir? 



Li-Xin Li 

Max-Planck-Institut fur Astrophysik, 85741 Garching, Germany 
e-mail: lxl@mpa-garching.mpg.de 

Abstract. One of the most important discoveries in the observation of gamma- ray bursts (GRBs) 
is that the total energy emitted by a GRB in y-rays has a very narrow distribution around 10 51 erg, 
which has led people to claim that GRBs are standard energy explosions. As people made the claim 
they have ignored the selection biases which must be important since GRB observations are strongly 
fluence or flux-limited. In this paper we show that, when the selection effects are considered, the 
intrinsic distribution of the GRB energy can be very broad. The number of faint GRBs has been sig- 
nificantly underestimated because of the fluence or flux limit. The bright part of the distribution has 
been affected by another important selection effect arising from the beaming of GRB jets, which is 
instrument-independent and caused by the fact that brighter GRBs tend to have smaller jet angles and 
hence smaller probabilities to be detected. Our finding indicates that GRBs are not a standard energy 
reservoir, and challenges the proposal that GRBs can be used as standard candles to probe cosmology. 
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1 Introduction 

A characteristic observed feature of cosmological gamma-ray bursts (GRBs) is that they emitted a 
huge amount of energy in y-rays in a very short time and their isotropic-equivalent y-ray energy 
(i.e., the total y-ray energy emitted by a GRB if the GRB radiates isotropically) spans a very large 
range — more than five orders of magnitude. The measured isotropic-equivalent energy of GRBs, 
E iso , ap pears to have a log -no rmal d istribution with a mean ~ 10 53 erg, and a dispersion ~ 0.9 in 
log £ 1S0 jAmatil200d 12001 11112001 . 

However, the re is evidence that GRBs are beamed dHarrison et alii 19991 : iKulkarni et alll 19991 : 
IStanek et al.ll 19991) . Assuming that a GRB radiates its energy into two oppositely directed jets, each 
having a half-opening angle 8 ia . The total solid angle spanned by the jets is then 4na>, where w = 
1 - cosset < 1. If the emission of a jet is distributed more or less uniformly on its cross-section, the 
total y-ray energy emitted by the GRB is approximately E y = a)E iso , smaller than £ iso by a beaming 
factor a>. 

One of the most important discoveries in GRB observations has been that the value of E y has 
a very narrow distribution with a mean ~ 10 51 erg comparable to ordinary supernovae, which has 
led people to claim that GRBs are a standard energy reservoir involving an approximately constant 
explosion energy (Frail et al, 2001; Piran et al. 2001; Berger, Kulkarni & Frail 2003; Bloom, Frail & 
Kulkarni 2003; Friedman & Bloom 2005 ). Theoretical models for interpreting the clustering GRB 
energy have also been proposed (see, e.g. JZhang & Me szaros 20021). 

It is well-known that observations of GRBs are strongly fluence or flux-limited, hence GRB 
samples seriously suffer from Malmquist-type selection biases dMalmquisill 1 920l ; iTeerikorpill 1 99l . 
That is, an observer will see an increase in the averaged luminosity or the total energy of GRBs with 
the distance, caused by the fact that less luminous or sub-energetic bursts at large distances will not 
be detected. Although this Malmquist bias for a flux-limited sample of astronomical objects looks 
obvious, sometimes people made serious mistakes in interpreting data by neglecting it. For instance, 
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with a study of n earby galaxies it had been incorrectly claimed that the Hubble constant increases 
with the distance Jde Vaucouleursll 1972k iTeerikorpill 1975k |Sandagelll994l) . 

Unfortunately, as people drew the conclusion on the distribution of the GRB energy and claimed 
that GRBs are standard energy explosions, they have treated the observed distribution as the intrinsic 
distribution and have neglected the selection biases that are very important for GRBs at cosmological 
distances. As a result, the number of faint GRBs has been significantly underestimated, since a GRB 
will not be detected if its flux or fluence falls below the detection limit. The bright part of the 
GRB energy distribution suffers from another important selection bias, which arises from the fact 
that brighter GRBs tend to have smaller jet opening angles and hence smaller probabilities to be 
detected. This beaming bias is independent of instruments and thus cannot be reduced by improving 
the sensitivity of detectors. 

The aim of this paper is to show that the influence of the selection biases from the fluence limit 
and the jet beaming is strong enough that the observed distribution of the GRB energy does not 
represent the intrinsic distribution at all. We present a simple model that explains nicely the observed 
distribution of the GRB energy, yet the burst energy reservoir in the model is not standard. Hence, 
the collimation-corrected energy of GRBs can have a very broad intrinsic distribution despite the fact 
that it is observed to cluster to a narrow distribution. Our results lead to the suggestion that GRBs 
are not a standard energy reservoir, contrary to the previous claim. 



2 Intrinsic versus Observed Energy Functions of GRBs 

To estimate the influence of selection biases on the observed distribution of the GRB energy, we 
assume that a GRB will be detected if one of its jets points toward the observer, and its fluence 
exceeds the limit F bolJim = 1.2 x 10~ 6 erg cm" 2 . Although this is an over-simplified approximation 
for the selection effect for GRBs, we will see that this simple selection effect already affects the 
observed distribution of the GRB energy strongly enough. In addition, it appears that the above 
fluence limit can reasonably represent the selec tion effect for GRBs with measured peak spectral 
energy and isotropic-equivalent energy dLil2007l) . 

The limit in fluence corresponds to a lower limit in the isotropic-equivalent energy of a de- 
tectable GRB at redshift z. -Bisojim = 4^-D^ om (l + z)-F bo ij im , where D com is the comoving distance to 
the burst. Here we assume a cosmology with £l m = 0.3, Da = 0.7, and a Hubble constant Hq = 70 
km s~' Mpc~'. 

For the GRB rate, as peop le often do we adopt the simplest as sumption that GRBs follow the 
cosmic star formation history dTotanilll997l ; iNataraian et"aUl2005h . Then, up to a normalization 
factor, the intrinsic distribution of GRB redshifts is given by 

j., n 2-sfr(z) dV com 

f(z) = — — , (1) 

1 +Z dz 

where S S fr(z) is the comoving star formation rate, and V com is the comoving volum e. 
We adopt a star formation rate dHopkins & Beacoml2 006: Le & Dermerll20071) 

v ( \ 1 +az rr\ 

Wz) = TTW (2) 

The parameters a, b and c are not well constrained. However, a model with a = 8, b = 3 and c = 1.3 
fits the observed distribution of GRB redshifts reasonably well dLe & Dermerll 20071) . Hence, we fix 
a, b and c to these values. 

For simplicity, we assume that except the number density, the pro perty of G RBs does not evolve 
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with the cosmic redshift, although this might not be true in reality B l2007h . Then, the intrinsic 
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distribution function of z, E- 1B0 and y = log (tan 0j el ) must have a form 

P(z,£W>,y) = f(z)(p iS o(E iso )^(E iso ,y) , (3) 

where we assume that \jj(E lso ,y) is normalized with respect to y: J {j/(Ei m ,y)dy = 1. We choose 
E is0 rather than E y as an independent variable, since in practice £ iso is easier to measure than E y 
although the latter might be more fundamental. 

For a GRB with a beaming factor <±>, the probability for it to be detected by an observer is a), 
without consideration of the fluence limit. The observed distribution of £ iso is then 

(f>im (£iso) = 0iso (£iso) <W> (£ is0 ) B (£ iso ) , (4) 

where 

<<jj> (E iso ) = J u)ip(E lm ,y) dy (5) 
is the averaged beaming factor, and the function 

S(E iso ) = f(z)dz (6) 



reflects the selection effect from the fluence limit. 

In equation ((6), for a given E is0 , the value of z ]im = Zy im (E iso ) is solved from the equation 
E[ so = £i SO ,ii m , or just given by the maximum redshift of GRBs if E-, m > £i so- ii m at z = Zmax (assuming 
that the distribution of GRB redshifts is cut off at z = Zmax). 

The intrinsic distribution function of the collimation-corrected energy E y = <±>E iso , derived from 
the distribution in equation <|3j, is 

<p y (£ r ) = J or^iso (<^~ l E y ) ^(w^E^y) dy . (7) 



The observed distribution of E y is then 



y (£ r ) = J 0, so (w 1 £ y )i/'(w 'E y ,yjE^La l E y ^dy. 



(8) 



It i s observed that the jet opening angle of GRBs is anti-correlate d to the isotropic-equivalent 
energy jFrail et alJlioOlklBloom et alj|2003l : iFriedman & BloomllloO^) . We define x = log£ iso and 
assume that if/(x,y) has a Gaussian form 



i/r(x,y) = J- exp 



(y - mx - p) 2 



2crj 



(9) 



where m, p, and cr y are constants. 

For a given x, the normalized observed distribution of y is 

W)L = -™ = — . . . , • (10) 

J _ oo co(y)if/(x,y)dy <w}(x) 

Because of normalization, the selection effect from the fluence limit [the function E (£i so ) defined by 
eq. |6) is canceled out in equation d I ( >b so it does not influence the observed distribution of the jet 
opening angle for a given E iso . However, the selection effect from beaming [i.e., the function a>(y)] 
is retained. 
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Fi gure 1 ! The observed distribution of the GRB energy is shaped by the fiuence limit of the 
detector and the distribution of the jet opening angle. The solid curve shows the fluence-selection 
function defined by equation $6^ normalized by S = C f(z)dz, where we have set z max = 10. The 
dashed curve shows the averaged jet beaming factor, defined by equation l(5j- The detection criteria 
are defined as follow: a GRB is detected if (1) one of its jets points toward the observer; and (2) its 
observed fiuence exceeds the limit F hol j im = 1.2x 1CT 6 erg cirr 2 . The fluence-selection function leads 
to a quasi-exponential cut-off to the GRB energy function at the low-energy end. The jet beaming 
factor reduces the number of detected GRBs of high energy, caused by the fact that brighter GRBs 
tend to have smaller jet angles. 

A maximum-like lihood fit of equation \10$ [with \ff{x,y) given by eq. |9) to the 23 GRBs with 
available E- KO and 0j e t dFriedman & Bloom 2005) leads to m = -0.216, p = -0.825 and <x v = 0.148 
(£ iso in 10 52 erg). 

In Fig.Q] we show the fluence-selection function defined by equation ® (with z max = 10) and the 
averaged lj defined by equation l[5j- The fluence-selection effect affects the observed distribution of 
faint GRBs, while the beaming effect affects the observed distribution of bright GRBs dramatically. 
For a given intrinsic distribution of the GRB energy, the combination of these two effects determines 
the shape of the distribution observed by an observer (if other selection effects are neglected). 

Finally, we assume that the intrinsic function of the isotropic-equivalent energy of GRBs is a 
power law with an exponential cut-off 

(fee (E h0 ) = Ef so exp (-E iso IE J , (11) 

where a and E* are constant parameters to be determined from observational data. 
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Figure 2: The dotted line histogram is the observed distribution of E- lso for 48 long-duration GRBs, 
with the number of GRBs in each bin indicated by a dark point with Poisson error bars. The dashed 
curve shows the intrinsic distribution of log E im defined by a power law with an exponential cut-off 
in equation i ll It . with a = -0.733, and = 3.21 x 10 54 erg. The solid curve is the observed 
distribution of log£j so derived from equation (|4), which well fits the observed data, with^ = 0.49. 
The fluence limit of the detector is 1.2 x 10~ 6 erg cm 4 . The maximum GRB redshift is set to be 
Zmax = 6, in accordance with the redshift distribution of the 48 GRBs. The intrinsic distribution is 
normalized so that the area under the dashed curve is the same as that under the solid curve. 

3 Results 

Fitting the observed distribution of E iso for 48 long-duration GRBs d AmatfeOOd 120071 : IlH|2007|) by 
the function in equation l(4j, we get a = -0.733 and E+ = 3.21 x 10 54 erg (Fig.|2j. The reduced 
chi-square of the fit is^J = 0.49, indicating a very good fit. This result clearly shows the fact that the 
intrinsic distribution of GRB energy is very different from the observed distribution, because of the 
strong selection effects from beaming and fluence-limit. The observed distribution has a Gaussian 
shape, but the intrinsic distribution is consistent with a power-law with an exponential cut-offQ The 
derived intrinsic distribution indicates the existence of a large amount of faint GRBs, which have not 
been detected because of the fluence-limit. In addition, the intrinsic distribution Ei so (f>i so (E lso ) peaks 
at E lm ~ 10 54 erg, an order of magnitude larger than the value ~ 10 53 erg directly inferred from the 
observed distribution. 

With a and fixed at the above values, we then fit the observed distribution of E y for 23 GRBs 

'We tried to fit the observed distribution of E iso with a Gaussian intrinsic distribution of log E iso , 
but we obtained an unrealistically large mean of E isB ~ 10 60 erg. 
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Figure 3: Distribution of the collimation-corrected energy of GRBs, E r , derived from the intrinsic 
distribution function of the isotropic-equivalent energy in equation Jilt , with a = -0.733 and = 
3.21 x 10 54 erg (the dashed curve in Fig.|2j. The solid curve is the observed distribution of logE,, 
(eq.[8]l. The dashed curve is the intrinsic distribution of log E y (eq.[7]l. The dotted line histogram is 
the distribution of the measured E y for 23 long-duration GRBs, with the number of GRBs in each 
bin indicated by a dark point with Poisson error bars. By varying the normalization, the derived 
distribution (the solid curve) fits the observation perfectly, with^ = 0.26. The intrinsic distribution 
is normalized so that the area under the dashed curve is the same as that under the solid curve. 



dFriedman & Bloomf 2005 ) by equation {§}> varying only the normalization. The best fit is shown in 
Fig- El with;^ = 0.26. Although this is not an independent fit given the fact that £j so and #j e , are 
anti-correlated which has been adopted by our model, the goodness of the fit in is still impressive, 
confirming that the relation between E[ so and 9j a assumed in equation ((9j( is a good approximation. 

Fig. [3] shows the dramatic difference between the intrinsic distribution and the observed distri- 
bution of log£ y . The observed distribution (the solid curve) has an exponential decay toward the 
faint burst end, but the intrinsic distribution (the dashed curve) decays toward the faint end by a 
power law: E r 

^iso^iso (^iso) ' 



, (£ y ) oc f^- 45 [slightly faster than the decay of the intrinsic distribution of log £j s 
£°- 27 l. The existence of a large amount of undetected faint GRBs broadens the 

ISO J o 



intrinsic distribution of log E y significantly. 

Here, by 'faint GRBs' we refer to those bursts with an £"; s 



or £ v that is smaller than the E,.„ or 



E y at the maximum of the distribution of log£"; so and log£ y . Although the detection of highly sub- 
luminous and sub-energetic nearby GRBs 98042 5, 031203 and 060218 has led peop l e to propose that 
there exists a unique population of faint GRBs i Cobb et aTll 20061 : IPian et al]l2006l : ISoderberg et alj 
120061 : iGuetta & Delia Valid [20071 : iLiang et alj|2007l) , the results in this paper do not rely on the 
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Figure 4: This figure shows how the observed distribution of GRB energy sensitively depends on 
the selection biases from the detector fiuence limit and the beaming of GRB jets. Upper panel: The 
intrinsic (solid curve) and the observed (dashed and dotted curves) distribution function of the GRB 
isotropic-equivalent energy, £j so . The solid curve is a plot of £; so 0i so (£"i S oX defined by equation dl lb 
with a = -0.733 and = 3.21 x 10 54 erg. The dashed curves are plots of £j so 0j so (E- lso ), calculated 
by equation l[4j with z max = 10 and F bol Jim = 10~ 6 , 10~ 7 and 10~ 8 erg cirr 2 respectively (from right 
to left, counted by the peak). The dotted curve shows ,Ej so 0j so (E lso ) in the limiting case of Fboijim = 
(i.e., GRBs can be detected down to any small value of fiuence). All distributions are normalized 
so that the integral over E iso is unity. Lower panel: Similar to the upper panel but for the intrinsic 
(solid curve) and observed (dashed and dotted curves) distribution of the collimation-corrected GRB 
energy, E y = u)E lso , calculated by equations and l[8j with different values of Fb iji m as in the upper 
panel. All distributions are normalized so that the integral over E y is unity. 
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existence of this unique population of faint GRBs. The extension of the derived energy function of 
normal GRBs to the low-energy end already significantly broadens the GRB energy function. 

The effect of the selection biases is more clearly illustrated in Fig. [4] which shows the depen- 
dence of the shape of Ei S0 ^ so CEiso) (upper panel) and E y <j> y (£ y ) (lower panel) on the beaming effect, 
and how the shape changes with the fluence limit of the detector. If the fluence limit decreases by 
a factor 10 from 10~ 6 erg cirT 2 (or from 10~ 7 erg cm -2 ), the width of E lso (j>i so (£i s0 ) and E y tj> y {E-^ 
increases by ~ 0.8 in log£i so and ~ 0.5 in log£ r , respectively. Correspondingly, the value of E 1S0 at 
the maximum of E iso (f> iso (E is0 ) decreases by a factor ~ 10 08 , and the value of E y at the maximum of 
E y cf> y (_E y ) decreases by a factor ~ 10 ' 5 . 

If we had an ideal detector that can detect an arbitrarily faint GRB, we would see a distribution 
of energy given by the dotted lines in Fig. [4] which is dramatically different from the intrinsic dis- 
tribution. For example, the observed distribution of E y (Ei so ) would peak at ~ 10 48 erg, rather than 
~ 10 51 erg (~ 10 54 erg) as indicated by the intrinsic distribution. 



4 Conclusions and Discussion 

The observation that the total energy emitted in y-rays by long-duration GRBs clusters around 10 51 
erg dFrail et alj|2001 ), wh ich has been considered as the most intriguing finding in GRB research 



( Zhang & Meszaroj 



20041) . is only a superficial result since the strong selection biases from the 



detector selection effect and the beaming of GRBs have been ignored. The previ ous claim that the 
energy output of the central engine of long-duration GRBs has a universal value dFrail et aT1l200ll ; 
IPiran et alfcOOll) . which was derived from the above superficial result, is likely to be incorrect since 
the observed narrow distribution of E y is consistent with a broad intrinsic distribution of E y . 

In fact, our results show that for both log is iso (Fig.0 and log E y (Fig. [3}, the distribution on the 
left-hand side to the maximum is well modeled by the cut-off from the fluence limit (the solid curve 
in Fig.Q}. It would be surprising that the intrinsic distribution happens to have a low-energy cut-off 
that is coincident with the fluence limit cut-off. 

The influence of the flux or fluence limit of detectors on the observation of GRBs is well- 
known and has been taken into account either thoroughly or partly in many GRB works, e.g. in 
deriving the luminosity function of GRBs (Schmidt 1999, 2001; Firmani et al. 2004; Guetta, Piran & 
Waxman 2005; Liang et al. 2007). However, the influence has sometimes been ignored or seriously 
underestimated. The claim that the collimation-corrected energy of GRBs has a narrow distribution 
and hence GRBs are a standard energy reservoir is an example where the selection effects have been 
ignored and wrong physical conclusions have been drawn. 

Although it is generally conceived that the jet opening angle is anti-correlated to the GRB en- 
ergy, in the study on the luminosity functio n of GRBs the effe ct of jet beaming has o ften not been 
properly taken into account. For example, in lGuetta et al. I J2005h and lLiang etaLU2007h . an isotropic 
luminosity function was derived by comparing the model prediction with the observed flux or the lu- 
minosity distribution without a consideration of beaming, then the derived luminosity function was 
used to calculate a weighted and averaged beaming factor. As we can see from equation ((4), the 
isotropic luminosity function derived by them should be the product of the intrinsic isotropic lumi- 
nosity function and the averaged beaming factor as a function luminosity, not the intrinsic isotropic 
luminosity itself. 

Our res ults also challenge the proposal that GRBs can be used as s tandard candles to probe 
cosmology dBloom et alj|2003l ; iFriedman & Bloonill2005l ; ISchaeferll2007l and references therein). 
Although the constancy of the GRB energy is not a necessary condition for GRBs to be standard 
candles because of the identification of several good correlations among GRB observables, the exis- 
tence of a large amount of faint bursts that have not been observed might significantly increase the 
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scatter in those correlations. In addition, the recent work of lButler et al.l d2007l) indicates that some 
of those relations arise from partial correlation with the detector threshold and hence are unrelated to 
the physical properties of GRBs. 

Finally, a prediction of this work that can be tested with future observations is that as the sen- 
sitivity of GRB detectors increases the observed distribution of E y broadens towards the low-energy 
end. 
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